Sonoluminescence as a QED vacuum effect: Probing Schwinger's proposal 
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Several years ago Schwinger proposed a physical mechanism for sonoluminescence in terms of 
photon production due to changes in the properties of the quantum-electrodynamic (QED) vac- 
uum arising from a collapsing dielectric bubble. This mechanism can be re-phrased in terms of the 
| Casimir effect and has recently been the subject of considerable controversy. The present paper 

probes Schwinger's suggestion in detail: Using the sudden approximation we calculate Bogolubov 
£SJ ' coefficients relating the QED vacuum in the presence of the expanded bubble to that in the pres- 

^ , ence of the collapsed bubble. In this way we derive an estimate for the spectrum and total energy 

emitted. We verify that in the sudden approximation there is an efficient production of photons, 
and further that the main contribution to this dynamic Casimir effect comes from a volume term, 
as per Schwinger's original calculation. However, we also demonstrate that the timescales required 
to implement Schwinger's original suggestion are not physically relevant to sonoluminescence. Al- 
though Schwinger was correct in his assertion that changes in the zero-point energy lead to photon 
production, nevertheless his original model is not appropriate for sonoluminescence. In other work 
we have developed a variant of Schwinger's model that is compatible with the physically required 
timescales. 
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I. INTRODUCTION 



In this paper we shall concentrate on Schwinger's original proposal regarding sonoluminescence Jl|-0, that of 
photon production associated with changes in the QED vacuum state. His idea was to explain the sonoluminescence 
phenomenon, which consists in light emission by a sound-driven gas bubble in fluid ||, in the framework of the 
so-called dynamical Casimir effect. Our first aim is to verify, in a dynamic framework, that a sudden change in 
bubble size will cause efficient photon production, thereby indicating the possibility of an a priori interesting role 
for the dynamic Casimir effect in this condensed matter context. While we demonstrate that the key features of 
Schwinger's calculations are correct, this study also demonstrates that for other reasons (to do with the observed 
timescale of the phenomenon) the original approach of Schwinger is not physically relevant to sonoluminescence. In 
related work [p|-p^[ we have developed a different implementation of Schwinger's ideas regarding sonoluminescence 
that is compatible with the physically observed timescales. 

The idea of a "Casimir route" to sonoluminescence was developed by Schwinger in a series of papers One 
key issue in Schwinger's model is simply that of calculating static Casimir energies for dielectric spheres — and there 
is already considerable disagreement on this issue. A second and in some ways more critical question is the extent to 
which a change in static Casimir energies might be converted to real photons during the collapse of the bubble — it 
is this issue that we shall address in this paper. We estimate the spectrum of the emitted photons by calculating an 
appropriate Bogolubov coefficient relating the two states of the QED vacuum. 

Another model associating sonoluminescence with QED vacuum changes is the variant of Schwinger's proposal due 
to Eberlein [O [L5|. In contrast to Schwinger's quasi-static approach, Eberlein's model is truly dynamical but uses 
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a radically different physical approximation — the adiabatic approximation. The two models should not be confused. 
See @ for a deeper discussion of Eberlein's approach to sonoluminescence. 

Considerable confusion has been caused by Schwinger's choice of the phrase "dynamical Casimir effect" to describe 
his model. In fact, the original model is not dynamical and is better described as quasi-static as the heart of the 
model lies in comparing two static Casimir energy calculations: that for an expanded bubble with that for a collapsed 
bubble. In a series of papers Schwinger showed that the dominant bulk contribution to the Casimir energy of a 
bubble (of dielectric constant (inside) hi a dielectric background (of dielectric constant (outside) is 

4tt 3 f K Ank 2 dk i f 1 1 

^cavity — + 2-— il / ~~ 77j T^ - ~TlCK 
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= + -LncR 3 K*(—^-—L=)+---. (1) 
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There are additional sub-dominant finite volume effects [|16|-|18[. The quantity if is a high-wavenumber cutoff that 
characterizes the wavenumber at which the dielectric constants drop to their vacuum values. This result can also be 
rephrased in the clearer and more general form as J16[- 



/d 3 k 1 
(27r)3 2^ \ L ° inslde ^ ~ U outside{k)\ +'•• (2) 

where it is evident that the Casimir energy can be interpreted as a difference in zero point energies due to the different 
dispersion relations inside and outside the bubble. 

In contrast, Milton [fill, and Milton and Ng [p0|,f2l| strongly criticize Schwinger's result. Using what is to our minds 



a physically dubious renormalization argument leads them 21 to discard both the volume and even the surface term 



and to claim that the Casimir energy for any dielectric bubble is of order E he/ R. 

In |l6|-|ilf| an extensive discussion on these topics is found. Therein it is emphasized that one has to compare two 
different geometrical configurations, and different quantum states, of the same spacetime regions. In a situation like 
that of Schwinger's model one has to subtract from the zero point energy (ZPE) for a vacuum bubble in water the 
ZPE for water filling all space. It is clear that in this case the bulk term is physical and must be taken into account. 
In the situation pertinent to sonoluminescence, the total volume occupied by the gas is not constant (the gas is truly 
compressed) , and it is far too naive to simply view the ingoing water as flowing coherently from infinity (leaving voids 
filled with air or vacuum somewhere in the apparatus). Since the density of water is approximately but not exactly 
constant, the influx of water will instead generate an outgoing density wave which will be rapidly damped by the 
viscosity of the fluid. The few phonons generated in this way are surely negligible. Surface terms are also present, and 
eventually other higher order correction terms, but they prove to not be dominant for sufficiently large cavities [18|. 



II. BOGOLUBOV COEFFICIENTS 



As a first approach to the problem we study in detail the basic mechanism of particle creation, and test the 
consistency of the Casimir energy proposals previously described. With this aim in mind we consider the change in 
the QED vacuum associated with the collapse of the bubble, by keeping fixed the refractive index both of the gas and 
of the water. For the sake of simplicity we take, as Schwinger did, only the electric part of QED, reducing the problem 
to a scalar electrodynamics. Moreover, at this stage of development, we are not concerned with the dynamics of the 
bubble surface. In analogy with the subtraction procedure of the quasi-static calculations of Schwinger [fij-fij, and of 
Carlson et al. [p~6|— p~8|| , we shall consider two different configurations of space. An "in" configuration with a bubble of 
dielectric constant (-inside (typically vacuum) in a medium of dielectric constant (outside-, and an "out" one in which one 
has just the latter medium (dielectric constant (outside) filling all space. Strictly speaking we should compare a large 
bubble having radius i? m ax with a small bubble of radius i? m in- We are approximating the small bubble by zero volume 
on the grounds that the small bubble that is relevant to sonoluminescence is at least a million times smaller than 
the large bubble at the expansion maximum. Keeping i? m i n finite significantly complicates the calculation but does 
not give much more physical information. The above "in" and "out" configurations will correspond to two different 
bases for the quantization of the field. The two bases will be related by Bogolubov coefficients in the usual way. Once 
we determine these coefficients we easily get the number of created particles per mode and from this the spectrum. 
This tacitly makes the "sudden approximation" : Changes in the refractive index are assumed to be non-adiabatic, 
see [p|-|i"l[ for more discussion. We shall also make a consistency check by a direct confrontation between the change 
in Casimir energy and the direct sum, E = WfcTifc of the energies of the created photons. The former energy (the 
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total energy of the particles that can be produced by the collapse) must necessarily equal the Casimir energy of the 
bubble in the "in" state since in the current simplified model there is no external source of energy (like the driving 
sound in the true dynamical effect). For this reason we expect to be able to give a definitive answer on the nature 
(dependence on the bubble radius and on the cut-off) of the static Casimir energy. Of course it is evident that such a 
model cannot be considered a fully satisfactory model for sonoluminescence. In fact this model completely ignores the 
details of the dynamics and moreover, by considering just one cycle, implies impossibility of testing for the possible 
presence of any parametric resonances. We thus consider the present calculation as a toy model in which some basic 
features of the Casimir approach to sonoluminescence are investigated: it provides a test of the nature and quantity 
of the particles produced by a collapsing dielectric bubble in the sudden approximation. 



A. Formal calculation 



Let us consider the equations of the electric fields (Schwinger framework) in spherical coordinates and with a time 
independent dielectric constant (we temporarily set c = 1 for ease of notation, and shall reintroduce appropriate 
factors of the speed of light when needed for clarity) 

ed (doE) - V 2 E = 0. (3) 

We look for solutions of the form 

E = f>(r,t)Y lm (^. (4) 



Then one finds 



e(d 2 $) - (d 2 $) + h(l + 1)$ = 0. (5) 



r 

For both the "in" and "out" solution the field equation in r is given by: 

edfo - 9 2 $ + + I)$ = 0. (6) 

In both asymptotic regimes (past and future) one has a static situation (either a bubble in the dielectric, or just the 
dielectric) so one can in this limit factorize the time and radius dependence of the modes: $(r, t) = e ""'/(/•). One 
gets 

/" + - ^1(1 + 1)) / = 0. (7) 

This is a well known differential equation. To handle it more easily in a standard way we can cast it as an eigenvalues 
problem 

/" - + 1)) / = -A 2 /, (8) 



where A 2 = ecu 2 . With the change of variables / = r 1 / 2 G we get 



G" + -G + ( A 2 - ^ ) G = 0. (9) 



r 

This is the standard Bessel equation. It admits as solutions the first type Bessel and Neumann functions, J v {Xr) and 
N v (\r), with v = I + 1/2. Remember that for those solutions which have to be well-defined at the origin, r = 0, 
regularity implies the absence of the Neumann functions. For the "in" basis we have to take into account that the 
dielectric constant changes at the bubble radius (R). In fact we have 

{^inside = n l as = dielectric constant of air-gas mixture if r < R, . . 

^-outside ^liau id " dielectric constant of ambient liquid (typically water) if r > R. 

Typically one simplifies calculations by using the fact that the dielectric constant of air is approximately equal 1 
at standard temperature and pressure (STP), and then dealing only with the dielectric constant of water (nu qu id = 
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yjtoutside ~ 1.3). We find it convenient to explicitly keep track of n gas and nu qu id in the formalism we develop. 
Defining the in and out frequencies, u)i n and u> ut respectively, one has 

r in (n n, .1,1. r\ — I A v Jv\n g a,sWi n r) if r < R, 

u i/ V l gasi 'l'hqmdi UJ ini ' ) — "i d 7 / „\ , AT(n,. .,r,|. t-"| if r /? V ^ 

The ^4 W , and C„ coefficients are determined by matching conditions in i? 

A v J v (rig as LOi n R) = BuJv{ n liquid^inR) ^ C u N v {niiq U i ( iLOi n R), ("121 

A v JJ{n gas Lo in R) — B u J u '(ni iqui dUJinR) + C u N u '(ni iquic iLUinR). 
The "out" basis is easily obtained solving the same equation but for a space filled with a homogeneous dielectric, 

G1 ut (nii qui d,uJout,r) = J v {nu qu idLQ out r). (13) 
To check that the "out" basis is properly normalized we use the scalar product, defined as usual by 

(&,&) = -i / 0i do<t>* 2 d 3 x. (14) 
is 

There are subtleties in the definition of scalar product which are dealt with more fully in The naive scalar 

product adopted here is missing a dependence on the refractive indices of the gas and the surrounding water. Given 
the fact that in the present framework both of these are approximately equal to one, the product adopted here is 
good enough for a qualitative discussion. Consider now the scalar product of a eigenfunction with itself, one expects 
to obtain a normalization condition which can be written as 

(«*)*> (KutT) = 5 ij - (15) 
Inserting the explicit form of the $ functions we get 

((<0*> (Ct)*) = (A + A') / rdrMXr)U\'r)e^-^ (16) 



= (A + y) (?(A 7 A,) e t(A - A ' ) , (17) 
A 

where we have used the Hankel Integral Formula 

i 

rdrJu{\r)J v {\'r) = 8{X - A')/A. (18) 
The Bogolubov coefficients are defined as 
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/% = +(£° ut ,£™*). (20) 

We are mainly interested in the coefficient /3, since |/3| 2 is linked to the total number of particles created. By a direct 
substitution it is easy to find the expression: 

f 3 = - l J Q f*«*(r,i) Yi m (Q) ^ do f$i„(r,f) *W(0) ^ r 2 drdf>, (21) 

= (u) in - uj out ) e l ( Uout+u;in)t 8w S mm > / G° ut (nu quid ,uj ou t,r) G" l (n gQS , nu quid , uj in , r) rdr. (22) 



To compute the integral one needs some ingenuity. Let us write the equations of motion for two different values of 
the eigenvalues, A and /i: 



G A + -G x +^--(1 + -) 2 ) G\ = 0, (23) 

G; + iG; + f At 2 -ia + i) 2 )G^o. 124, 
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If we multiply the first by G M and the second by G\ we get 



G\G U H — G X G U + [ A 



G^G\ + —G^G\ 



2' 
2 J 



G\G^ — 0, 
G^Gx = 0. 



Subtracting the second from the first we then obtain 

1 



GxGu — G. G\ 



G A G /i - G.Gx ) + (A - ^)G X G„ = 0. 



The second term on the left hand side is a pseudo-Wronskian determinant 

W X)1 (r)=G' x (r)G^r)-G'^r)G x (r), 



(25) 
(26) 

(27) 
(28) 



and the first term is its total derivative dWx^/dr. It's a pseudo-Wronskian, not a true Wronskian, since the two 
functions G x and G M correspond to different eigenvalues and so solve different differential equations. The derivatives 
are all with respect to the variable r. Using this definition we can cast the integral over r of the product of two given 
solutions into a simple form. Generically: 



r-b r-h r-b 

(/i 2 - A 2 ) / rdr G A G M = / rdr dW x ^ + dr W X r 

J a J a J a 



That is 



So the final result is 



/' 

■J a 



rdr G A G M = 



1 



(M 2 - A 2 ) 



Wx,r 



f 

J a 



r dr GaG m = 



1 



- A 2 ) 

This expression can be applied in our specific case equation (| 
r dr G° ut (nu qmd ,uJ,r) G™{n gas ,ni iquid ,uj,r) 



drWx, 



(Wx^r) . 

a 

2), wc obtain: 



dr WV 



(29) 



(30) 



(31) 



r dr G° v ut (n Hquid LJ out r)G™(n gas u in r) + / r dr G° u ut (n Uquid uJ utr)G" l (nu qmd aj in r) 

JR 

{rW[G° ut (n Hquid uj out r), G™(n gas uw)]} Q {rW[G° ut (n Hquid uj out r), G'J 1 {n Uquid uj in r)}}™ 



{^liquid ^out)^ {j^gas ^in)^ 



illiquid bJout) illiquid ^inj 



R 



W[G v ut {n Uquld LJoutr),G l J l (n gas uj m r)] R _ W[G° ut (n Uqmd cj out r),G l J l (ni iqmd uj m r)] R+ 



illiquid &out) (j^gas ^in) 



illiquid ^out) {j^liquid <-^m) 



(32) 

(33) 
(34) 



where we have used the fact that the above forms are well behaved (and equal to 0) for r — 0, and r = 00 by 
construction. (Here and henceforth we shall automatically give the same I value to the "in" and "out" solutions by 
using the fact that equation (p2|) contains a Kronecker delta in I and I'.) 

Finally the two pseudo-Wronskians so found can be shown to be equal (by the junction condition (|l2|)). In fact one 
can easily check that: 



A u W[J u (n u 

quid ^out 

r), J v (n gas u>i n r)] R = B L ,W[J„( 

^liquid ^out^ r )]R 
+ C V W[J„{ ^liquid ^out 

r),N v (n gas Ui n r)]R. 



(35) 



This equality allows to rewrite integral in equation (E2J) in a more compact form 
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Jo 



r dr G° ut {n Hquid , ui, r) G" l (n gas ,n Hquid , cj, r) 
A 



quid ^out ) 2 - {ngas OJin ) 2 {niiquid ^out ) 2 ~ ^liquid UJzn ) 5 



W[J v (nii qu id uJoutr), J v {n g as uJ in r)} R (36) 



n Kquid- n 9 as\ A v Rwj n W[J v {nu qui d UJoutr), J v {n gas u) in r)]i 



Il- 



liquid 



\uJ 2 out - Ujf n ] [{nuquid L0 O ut) 2 - {n gas 0J in )' 



(37) 



Inserting this expression into equation ( p2| ) we get 

{^>out ~ Ui n ) 



= 



n liquid U gas 



Sw 5, 



liquid 



mm' 2 



RA u 2 n W[J v {n Uquld ujputr), J v {n gas uj m r)] R j(u out+Uin ) t 

[{niiquid U ou t) 2 — {n gas Wi„) 2 ] 



(38) 



We are mainly interested in the square of this coefficient summed over / and m. It is in fact this quantity that is 
linked to the spectrum of the "out" particles present in the "in" vacuum, and it is this quantity that is related to the 
total energy emitted. Including all appropriate dimensional factors (c, h) we have 



dN(u ou t) 
du> m ,,t 



\(3{u) in ,uj ut)\ 2 duj 

i 



and 



E = h 



dN{uj out ) 

dOJnut 



^outduJout- 



(39) 



(40) 



Hence we shall deal with the computation of: 

\(3{uj in ,U! ou t)\ 2 = ££ [Am,l'm'( w m!W otl( )] 2 



(41) 



lm I'm' 

/ 2 2 v 

/ ^liquid ^gas 



Il- 



liquid 



^out + 



2 oo 



£(21 + 1)1^,1 



i=i 



/quia. 

r/c)]i 



illiquid ^out) {j^gas ^inj 



(42) 



This expression is too complex to allow an analytical resolution of the problem. Nevertheless we shall show that the 
terms appearing in it can be suitably approximated in such a way as to obtain a computable form that shall give us 
some information about the main predictions of this model. We shall first look at the large volume limit, which will 
allow us to compare this result to Schwinger's calculation, and then develop some numerical approximations suitable 
to estimating the predicted spectra for finite volume. 



B. Behaviour in the large R limit 

One of the main objectives of this calculation is to shed some light on the extent to which the change in static Casimir 
energy can be transformed into photons. In particular we expect that the total energy of the photons calculated from 
this Bogolubov approach would give approximately the same results as the static Casimir energy calculations such 
those of Schwinger, and of Carlson et al. j|,[l9]|l(| , since we have excluded any external forces. 

From equation ( p2[ ) it is easy to check that the general form of the squared Bogolubov coefficient is given by 

\P{x,y)\ 2 = ^f3 2 {x,y), (43) 
cr 

where 0^{x,y) is a dimensionless quantity and we introduce dimensionless variables x ~ nii qu id^outR/c and y ~ 
n gas uiinR/c. (The dimensions of (3 should always be those of time.) The spectrum is then given by 

dN{u out ) R ( f RK , .. a , \ .... 

= / \Po(x,v)rdy , (44) 



dugout c n ga s 



and the energy radiated is 



G 



n liquid U 9as & 



/ dx dyx\fj (x,y)\ 2 . (45) 
Jo Jo 



If R is very large (but finite in order to avoid infra-red divergences) then the "in" and the "out" modes can both be 
described by ordinary Bessel functions 

G m (n gasi w,r) = J„(n gas (j in r/c), (46) 

G out (n Hquid ,uj,r) = J v (ni iquid uj out r/c). (47) 

We can now compute the Bogolubov coefficient relating these states 

Pa = (48) 



= (u m J^out) el ( Uout+0Jin )t Su> J J v (n gas u in r/c)J u (nu quid uj out r/c) r dr 

(, , , , \A(ui out +u) in )t r r 5{n gas LOjn — Uliguid ^>out) 
= {Win —Uout)e 0u> mm > 



n nn * lu 

1 1 



(49) 
(50) 

e imn(n gae /n liquid + l)t §w 6mm/ S ( Ugas ^ _ nHqmd ^ 

^gas ^liquid J 

This result implies that 

\P(uJ in , UJout)\ 2 = y^y] \f3lml'm'(Uin, Wout)\ 2 (52) 
Im I'm' 

2 



= — — _£2£j ^(21 + 1) {6(n gas u in - n Uquld u out )} 2 (53) 

\ illiquid Tlgas J . 

= ( nUqmd — ^£ J ^(21 + 1) 6(0) S(n gas u) in - n Uquid u) out ) (54) 

\ ^liquid Ugas J ^ 

= (— — IVM ^(21 + 1) — 6(n gas u in - nu quid w ou t), (55) 

\ n liquid IT-gas J ^ Z7TC 

where we have invoked the standard scattering theory result 

{5 3 (k)¥ = j^S 3 (k), (56) 

specialized to the fact that we have a I-dimensional delta function (in frequency not momentum). The sum over 
angular momenta (which is formally infinite) can now be estimated as follows 

'max(^out) 

J2 ( 2 * + 1 ) = CzK«t)-i~CsK u t). ( 57 ) 
i=i 

The angular momentum cutoff is estimated by taking 

, , ^ (illiquid hUout/c) X R 

lmax(u) ~ ^ = n Uquid UJ out R/C = X. (58) 

So in the above we are justified in approximating 

J2( 2l + 1 )~ x2 - (59) 



By changing to the dimcnsionless variables (x, y) this finally gives 

2 

^liquid Tlgas J 2,7TC^ 



\P(x,y)\ 2 = ( nii(luid n9as ) A- x 2 6(x-y). (60) 



7 



We can now compute the spectrum and the total energy of the emitted photons 
dN(uj ou t) R f Tinquid — n g as N r '" x 



di^out 'Z'KC Tlgas \ Tlliquid Tlgas J J Q 

2 



x 8{x — y) dy 

J Jo 

p / _ \ 2 

/ ^liquid Tlgas \ 2 



27TC Tlgas \ illiquid Tlgas 
^liquid ( ^liquid ^#as \ ^^^out 



) ( K - n ^^«*/ c )' ( 61 ) 



illiquid Tlgas 

For future comparison purposes it is convenient to write the spectrum in dimensionless form as 



dN If nii qu i d -n qi 



dx 27T Tliiq U id Tlgas \ TlHq U id Tlgas 

^From this equation it is easy to get the total number of the created photons 

1 / Tlliquid Tlgas 



x 2 G(RK-x). (62) 



N : 



27T Tlliq U id Tlgas \ Tlliquid Tlgas J JQ 



/ x 2 e(RK~x)dx 
Jo 



1 / Tlliquid Tig 



2 

gas \ ,p7^\3 



67T Tlliquid Tlgas \ TlHq U id Tlgas 

and the total emitted energy 



{RKy . (63) 



TlC l Tl IT \ ^ f°° fRK 

( nT^n ) / ^ x ^ x ^ ~ ^ 

^ "-liquid n 9 as 11 \ "liquid Ugas / Jo J 

* / _ \ 2 rRK 

/ TlUquid Tl, 



"/""' "" "•/«» I f 



^ n liquid n aa-s R \ Tln qu id Tlgas J Jo 

TlC ( nuq ui d - Tlg as \ ^ {RK) 4 



^ 7rn Uquid n 9asR \ Tlliquid Tlgas ) 4 



2 



n Uguid Tlgas\ hcK ( RK )3 ( g 4) 



""liquid n gas \ Tlu qu id Tlgas 

Hence, feeding our results ( |60| ) into equations (|44|) and ([45|) for dN{ui)/dui and E 1 , we find a result which is in 
substantial agreement with the Schwinger (and Carlson et al.) results. We view this is definitive proof that indeed 
Schwinger was essentially correct: The main contribution to the Casimir energy which can be extracted from the 
collapse of a (large) dielectric bubble is a bulk effect. The total energy radiated in photons balances the change in 
the Casimir energy up to factors of order one which the present analysis is too crude to detect. (For infinite volume 
the whole calculation can be re-phrased in terms of plane waves to accurately fix the last few prefactors.) 

In Schwinger's original model he took n gas ~ 1, Tiuq U id ~ 1-3, R ~ i?max ~ 40 /im and K ps 2tt/(360 nm) @. Then 
KR « 698. Substitution of these numbers into equation ([!]) leads to an energy budget suitable for about three million 
emitted photons. 

By direct substitution in equation (^3|) it is easy to check that Schwinger's results can qualitatively be recovered 
also in our formalism: in our case we get about 0. 7 million photons for the same numbers of Schwinger and about 
1.5 million photons using the updated experimental figures i? m ax ~ 45 /jm and K « 271"/ (300 nm). 

The average energy per emitted photon is approximately^ 

3 3 

(E) = -he K/niiq la d = ~h w max ~ 3 eV. (65) 

It is important to stress that equation ([!]) and equation ( |64| ) are not identical (even if in the large R limit the 
leading term of Casimir energy of the "in" state and the total photon energy coincide). One can easy see that the 



1 The maximum photon energy is h aj max ~ 4 eV. 
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volume term we just found [equation (J64|)] is of second order in (n — 1) and not of first order like equation (Q). This 
is ultimately due to the fact that the interaction term responsible for converting the initial energy in photons is a 
pairwise squeezing operator Q . Equation (|||) demonstrates that any argument that attempts to deny the relevance 
of volume terms to sonoluminescence due to their dependence on (n — 1) has to be carefully reassessed. In fact what 
you measure when the refractive index in a given volume of space changes is not directly the static Casimir energy 
of the "in" state, but rather the fraction of this static Casimir energy that is converted into photons. We have just 
seen that once conversion efficiencies are taken into account, the volume dependence is conserved, but not the power 
in the difference of the refractive index. 

Indeed the dependence of \f3\ 2 on (n — l) 2 and the symmetry of the former under the interchange of "in" and "out" 
state also proves that it is the amount of change in the refractive index and not its "direction" (from "in" to "out" ) 
that governs particle production. This also implies that any argument using static Casimir energy balance over a full 
cycle has to be used very carefully. Actually the total change of the Casimir energy of the bubble over a cycle would 
be zero (if the final refractive index of the gas is again 1). Nevertheless in the dynamical calculation one gets photon 
production in both collapse as well expansion phases. (Although some destructive interferences between the photons 
produced in collapse and in expansion are conceivable, these will not be really effective in depleting photon production 
because of the substantial dynamical difference between the two phases and because of the, easy to check, fact that 
most of the photons created in the collapse will be far away from the emission zone by the time the expansion photons 
would be created.) This apparent paradox is easily solved by taking into account that the main source of energy is 
the sound field and that the amount of this energy actually converted in photons during each cycle is a very tiny 
amount of the total power. 

We now turn to the study of the predictions of the model in the case of finite radius. Unfortunately this cannot be 
done in an analytic way due to the wild behaviour of the pseudo-Wronskian of the Bessel function. Nevertheless some 
ingenuity and a detailed study of the different parts of the Bogolubov coefficient leads to reasonable approximations 
that allow a clear description of the spectrum of particle predicted by the model. 



C. The A factor 



The A\, Bp, and C M factors can be obtained by a two step calculation. First one must solve the system ( |l2| ) by 
expressing B and C as functions of A. Then one can fix A by requiring B 2 + C 2 = I , a condition which comes from 
the asymptotic behaviour of the Bessel functions. Following this procedure, and again suppressing factors of c for 
notational convenience, we get 



A u = 



W[J v (nu quid u> in r), N v (illiquid u in r)} 



y/W[J v (n gas u> in r), N v (nuquid w in r)] 2 + W[J v (n gas u> in r), J v (n Uquid uj in r)} 2 
),N„(n H 

quid. 



A 



W[J v (n H 

quid ^in 

r),N v (nu 

quid LOinT )\ 

W[J v (nu qu id ui in r), J v (n gas uj in r)] 



W[J V (niiq U id Lo in r), N v (n liqmd uj in r)] 



(66) 

(67) 
(68) 



We are mostly interested in the coefficient A v . This can be simplified by using a well known formula (Abramowitz- 
Stegun, page 360 formula 9.1.16) for the (true) Wronskian of Bessel functions of the first and second kind. 



W true [J v (z),N v (z)\ 



(69) 



In our case, taking into account that for our pseudo Wronskian the derivatives are with respect to r (not with respect 
to z), one gets for the numerator of A v : 



W[J v (m u}i n r),N u (ni iquid ^m^)]i? — ^liquid 



^(illiquid Wi n R) 



2 

ttR' 



Hence the A,, can be written as 



\A„ 



4/(7T 2 i? 2 ) 



W[J u (n gas 0Ji n r),N u (nuq U id w,„r)] 2 + W[J u (n gas 0Ji n r), J v (niiquid Wmr)]\ 
For uj — ► oo at / fixed the asymptotic behaviour is 



(70) 



(71) 
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1^ |2 ^ 2n gas nii qu id {72) 

n 2 gas + nlquid + (rfiquid ~ n las) C0s(2n gas 0J in - [l> + 1/2)tt) ' 

Numerical plots of \A„\ 2 show that it is an oscillating function which rapidly reaches this asymptotic form. 




FIG. 1. Plot of | A;| 2 for v = 3/2 and v = 5/2. (Here we define z = ujR/c and put n gas — 1, nu qu id = 1.3.) The function 
rapidly stabilizes to the asymptotic behavior. The dotted curve shows the behaviour of the asymptotic form for v = 3/2. 



We shall use this approximation to replace the A v factor with its mean value for large arguments: 

,.2ji 



1-4, 



1 

2^ 



dz- 



2^-f/as ^liquid 



las + rfiquid + ( n Lmd ~ n jas) C0S ( Z ) 



= 1. 



(73) 



That this approximation is adequate may be checked a posteriori by seeing that the Bogolubov coefficients are not 
noticeably affected. 



D. The Pseudo-Wronskian 



Use the simplified notation in which x = nu qu id<-0 O utR/c, y — 



L gas win 



R/c. In these dimensionless quantities, after 



including the approximation equation (73), and making explicit the dependence on R and c, equation (42) takes the 
form: 



I , N ,2 R 2 ( n liquid n gas) 

\/3{x,y)\ - 



c2 n liquid n gas \n gas X + Illiquid y 



y 



F{x,y). 



(74) 



Here F(x, y) is shorthand for the function 



Jv{x) J v {y) 
xJ' v {x) yJ'M 



(x 2 - y 2 ) 



2)2 



(75) 



where in this equation the primes now signify derivatives with respect to the full arguments (x or y). 

In order to proceed in our analysis we need now to perform the summation over angular momentum. Although 
the infinite sum is analytically intractable, there are two reasonable arguments (one physical and one mathematical) 
both leading to the conclusion that suitable truncations of this sum will be enough for our purposes. 

The first argument is a physical one and it is based on the maximum amount of angular momentum that an 
outgoing photon may have. Basically, if one supposes the photons to be produced inside or at most on the surface of 
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the bubble, the upper limit for the angular momentum will be the product of the bubble radius times the maximal 
"out" momentum. Then one gets: 



lma X {K) = -^—^ = RK. 



(76) 



For sonoluminescence K is of order 27r/(200nm). Deciding the appropriate value of R is more tricky. Since the 
sonoluminescence flash occurs at or near the moment of minimum radius one might wish to use R m in ~ 500nm in 
which case l ma x{K) ~ 15. Certainly for this choice of R keeping the first ten or so terms will be sufficient. More 
conservatively, one might wish to choose R to be of order R a mbient ~ 4.5^tm in which case l ma x(K) rs 135. Keeping 
this number of terms in the series is already very unwieldy. Finally, in Schwinger's original version of the model it 
is the change in Casimir energy during the collapse all the way from maximum radius that is relevant, so perhaps 
one should use R m ax ~ 45^m. In this case l m ax{K) ~ 1350, and explicit summation of the series is prohibitively 
difficult. To handle these problems we develop a semi-analytic approximation to the sum which is sufficient for making 
numerical estimates of the spectrum. 

This argument can be bolstered by considering the large order expansion (y — > oo at fixed a;) of the Bessel functions. 
In this limit one gets |23||: 



J„(x) 



(77) 



This can be used to obtain the asymptotic form of the pseudo-Wronskian appearing in equation d73) . 



W v (x,y) 



Ju{x) J v (y) 
x J' v (x) y J' v (y) 

J v {x) J„(y) 
x J„+i(x) y J u +i(y) 

(x 2 - y 2 ) ( xy 



2ir(v) 1 / 2 (v + l) 3 / 2 \v{v + l) 



(I) 



2v+l 



(78) 
(79) 
(80) 



where we have used the standard recursion relation for the Bessel functions J' v [z) = vJ v {z) — zJ u+ i(z). This indicates 
that the sum over v is convergent: the terms for which (xy/is 2 ) < 1 are suppressed. Since, depending on one's views as 
to the appropriate value of R, x and y are at most of order 15, 135, or 1350 we deduce that the maximal contribution 
to the sum comes from a limited number of terms. 

Analytically, it is easy to see that the function F(x, y) is finite along the diagonal and goes smoothly to zero for 
x,y — > 0. To proceed to an actual computation of the predicted spectrum we need to develop an semi-analytic 
approximate form for this function by considering separately the behaviour along the diagonal x — y — and in the 
transversal direction x + y = constant. 



E. Working along the diagonal 



To study in more detail the behaviour of such a function in this zone one can perform a Taylor expansion of F{x 1 y) 
around x = y. 



v W v {x,y) 

lim ~i r = 

x^y (x — y) 



lim 

x—>y 



Jy{x) J y {y) 
x J' v {x) y J' v (y) 



lim 

x— >y 



{x-y) 

Ju{x) J u (x) + (x-y)J'„(x) 

x J' v {x) x J' v {x) + {x- y)[J' v {x) + x J'J(x)] 



(x-y) 

J Ax) J'M 

x J' v {x) J' v {x)+x J'l{x) 
= J v {x)[J' v {x) + x J'J(x)] -xJl(xf. 

The derivatives can be eliminated by using the well known recursion relations. 



(81) 

(82) 

(83) 
(84) 
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r W v {x,y) 

lim — — = J v (x) 

x~>y [x - y) 



~{v 2 -x 2 y 




— X 




X 





-J v (x) - J v+ i{x) 



= 2vJ v {x)J v+1 {x) - x [Jl(x) + J 2 u+l (x)] . 
For sake of simplicity we shall use an equivalent form of equation ( |86| ) where lower order Bessel function appear 

W u (x,y) 



lim 

x^y [x — y) 



2vJ u (x)J v -i(x) - x [Jl(x) + Jl-xix)] . 



(85) 
(86) 

(87) 



This result shows that, as expected, each term of F(x,x) is finite along the diagonal and equal to zero at x = y = 0. 
Moreover 



D(x) = F(x,x) =5^(2J+1) 



z=i 



{(2Z + l)Ji + i/ 2 (x)J,_i /a (x) -x Jf +1/2 (x) + Jf_ 1/2 (x) } 



This sum can easily be checked to be convergent for fixed x. [Use equation (|77|).] With a little more work it can be 
shown that 



lim^D(x) = n _ 2 . 



1 

2T 2 



The truncated function obtained after summation over the first few terms (say the first ten or so terms) is a long 
and messy combination of trigonometric functions that can however be easily plotted and approximated in the range 
of interest. Due to numerical artifacts, the function is not controllable near the origin, fortunately we have analytic 
information in that region — the function is very near to zero in the range (0, 1) for both "out" and "in" frequencies, 
and can be approximated by zero without any undue influence on the numerical results. A semi-analytical study led 
us to the approximate form of D{x) 



A confrontation between the two curves in the range of interest is given in the figure below. 



(89) 




FIG. 2. Plot of the exact D(x) against its approximated form in the range 1 < x < 30 

F. The factorization approximation 

To numerically perform the integrals needed to do obtain the spectrum it is useful to note the approximate factor- 
ization property 
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fKs) , f (£±!.£±i) c (£_!). m 

That is: to a good approximation F(x,y) is given by its value along the nearest part of the diagonal, multiplied by 
a universal function of the distance away from the diagonal. A little experimental curve fitting is actually enough to 
show that to a good approximation 

F(x v) « D ( sin2(7r[x ~^ ]/4) (91) 

From the plot we show below it is easy to check that the function F(x, y) is quite well fitted by our approximation. We 
feel important to stress that this is approximation is based on numerical experimentation, and is not an analytically- 
driven approximation. (In the infinite volume case we know that F(x,y) — ► (constant) x 5{x — y), cf. equation fefj). 
The effect of finite volume is effectively to "smear out" the delta function. The combination sin (x)/(ttx 2 ) is one of 
the standard approximations to the delta function.) Our approximation is quite good everywhere except for values 
of x and y near the origin (less than 1) where the contribution of the function to the integral is very small. 




FIG. 3. Transverse fit: An orthogonal slice of F(x,y) intersecting the diagonal at (x,y) — (3,3). Here F(3 + z, 3 — z) is 
plotted in comparison with [sm 2 (TTz/2)]/(ivz/2) 2 . 



U . 03 
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FIG. 4. Transverse fit: An orthogonal slice of F(x,y) intersecting the diagonal at (x,y) — (5,5). Here F(5 + 2,5 — z) is 
plotted in comparison with [sm 2 (nz/2)]/(-Kz/2) 2 . 




G. The spectrum: numerical evaluation 

We have now transformed the function F(x, y) into an easy to handle product of two functions 

^/ „n 1 (x + y-2) 2 sin 2 (7r[a; - yl/4) 

We exhibit tridimensional graphs for both the exact (apart from the approximation of truncating the sum at a finite I) 
and approximate forms of the function F(x, y). We have chosen the case of R — 0.5^m (corresponding to y max = 2.5 
as previously explained). 




FIG. 6. Plot of the exact F(x,y) in the range < x < 5, < y < 2.5 
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FIG. 7. Plot of the approximated F(x,y) in the range < x < 5, < y < 2.5 



A numerical study of the error due to the replacement of F(x,y) with its approximated form equation (|92j), leads 
to an upper limit of 20% error in the total energy emitted. 

The dimensionless spectrum, based on equations (f42|) and (174), is 



dN 
d.v 



( n liquid n gas) 



n liquid n gas 



RK 



y 



D 



x + y\ sin 2 (7r[a; — y]/4) 



dy. 



/,, \n gas x + nn qu i d y J V 2 J ( w [ x - 2/]/ 4 ) 2 
As a consistency check, the infinite volume limit is equivalent to making the formal replacements 



sin (7T \x 



tf]/4) 



{Tt[x-y]/Af 



and 



D 



x + y 



1 

2^2' 



(93) 



(94) 



(95) 



Doing so, equation (|93j ) reduces to equation (62) up to an overall factor [4/7r] of order one. The correct dependence 
on refractive index and correct power-law behaviour for the spectrum are recovered, and the overall order one factor 
is merely a reflection of the crudity of the cutoff in angular momentum used in deriving (|62|). 

With this consistency check out of the way, it is now possible to perform the integral with respect to y to estimate 
the spectrum for finite volume. For defmiteness we set nuquid — 1-3 and n gas = 1.0, put K = 27r/(200nm), and pick 
R = 0.5/im (corresponding to y ma x = 15). We integrate from y = to y — 15 and plot the resulting spectrum from 
x = to x = 18. 
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FIG. 8. Spectrum obtained by the approximated Bogolubov coefficient for R = 0.5/im corresponding to y max = 15. We 
integrate from y — to y — 15 and plot the resulting spectrum from x — to x — 18. 

One can also ask what sort of result one would get if instead we pick a much larger value of R, say R = 5//m, 
corresponding to the bubble at equilibrium radius. In this case the approach towards the Schwinger (infinite volume 
result) result is much closer. We now have y m ax = 135. We integrate from y — to y = 135 and plot the resulting 
spectrum from x = to x = 140. For comparison we plot it together with equation (62) which is Schwinger's naive 
model (the re-scaled infinite volume limit). 




FIG. 9. Spectrum obtained by the approximated Bogolubov coefficient for R = 5fim corresponding to y ma x = 135. We 
integrate from y = to y = 135 and plot the resulting spectrum from x = to x — 140. 

The case corresponding to the bubble at maximum radius, R — 50/im, requires a range of integration too large for 
standard numerical plotting. In any case the graph will only be a replica of the previous one on a larger scale. 

III. DISCUSSION 

The lessons we have learned from this test calculation are: 

— (1) The model proves (in an indirect way) that the Casimir energy produced via the bubble collapse includes (in 
the large R limit) a term proportional to the volume (actually to the volume over which the refractive index changes). 
In the case of a truly dynamical model one expects that the energy of the photons so created will be provided by other 
sources of energy (e.g., the sound wave), nevertheless the presence of a volume contribution appears unavoidable. 

— (2) The present model is still unable to fully fit other experimental features of sonoluminescence. For example 
it provides maximal photon release at maximum expansion. Barber et al. B point out that in Schwinger's original 
model the main production of photons may be expected when the the rate of change of the volume is maximum, which 
is experimentally found to occur near the maximum radius. In contrast the emission of light in sonoluminescence is 
experimentally found to occur near the point of minimum radius, where the rate of change of area is maximum. All 
else being equal, this would seem to indicate a surface dependence and might be interpreted as a true weakness of the 
dynamical Casimir explanation of sonoluminescence. 

In fact we have shown elsewhere |^-[l^] that the situation is considerably more complex than might naively be 
thought. It is important to stress that what Schwinger proposed was clearly only a first estimate of the vacuum 
energy, which was in principle viable as the basis for a model, and not a fully dynamical model. Schwinger was fully 
aware of this in his papers. 

A fully dynamical calculation is required in order to deal with these issues, and the experimental data give remarkable 
suggestions about the plausible directions for theoretical developments within the framework of the dynamic Casimir 
effect. In particular, one of the key features of photon production by a space-dependent and time-dependent refractive 
index is that for a change occurring on a timescale r, the amount of photon production is exponentially suppressed 
by an amount exp(— lot). In |lC|| we have provided a specific model that exhibits this behaviour, and argued that the 
result is in fact generic. The importance for sonoluminescence is that the experimental spectrum is not exponentially 
suppressed at least out to the far ultraviolet. Therefore any mechanism of Casimir-induced photon production 
based on an adiabatic approximation is destined to failure: Since the exponential suppression is not visible out to 
lo rs 10 15 Hz, it follows that if sonoluminescence is to be attributed to photon production from a time-dependent 
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dielectric bubble (i.e., the dynamical Casimir effect), then the timescale for change in the refractive index must be of 
order of a femtosecond. Thus any Casimir-based model has to take into account that it is no longer the collapse from 
-Rmax to i? m in that is important. One has to divorce the change in refractive index from direct coupling to the bubble 
wall motion, and instead ask for a rapid change in the refractive index of the entrained gases as they are compressed 
down to their van der Waals hard core |To|JTl| ] . We stress that this conclusion, though it moves away from the original 
Schwinger proposal, is still firmly within the realm of the dynamic Casimir effect approach to sonoluminescence. The 
fact is that the present work shows clearly that a viable Casimir "route" to sonoluminescence cannot avoid a "fierce 
marriage" between QFT and features related to condensed matter physics. 



IV. CONCLUSIONS 



The present calculation unambiguously verifies that a sudden change in radius of a dielectric bubble causes a change 
in the Casimir energy that is, as predicted by Schwinger [fl]]l6|-|l8|] , converted into real photons with a phase space 
spectrum. As far as sonoluminescence is concerned, we have also explained why such a change must be sudden in 
order to fit the experimental data. This leads us to propose a somewhat different model of sonoluminescence based on 
the dynamical Casimir effect, a model focussed on the actual dynamics of the refractive index (as a function of space 
and time), and not just of the bubble boundary. (In Schwinger's original approach the refractive index changes only 
due to motion of the bubble wall.) In summary, provided the sudden approximation is valid, changes in the refractive 
index will lead to efficient conversion of zero point fluctuations into real photons. Trying to fit the details of the 
observed spectrum in sonoluminescence then becomes an issue of building a robust model of the refractive index of 
both the ambient water and the entrained gases as functions of frequency, density, and composition. Only after this 
prerequisite is satisfied will we be in a position to develop a more complex dynamical model endowed with adequate 
predictive power. 

In light of these observations we think that one can also derive a general conclusion about the long standing debate 
on the actual value of the static Casimir energy and its relevance to sonoluminescence: Sonoluminescence is not 
directly related to the static Casimir effect. The static Casimir energy is at best capable of giving a crude estimate 
for the energy budget in sonoluminescence. 
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